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I. INTRODUCTION 

Density Functional Theory (DFT) [l|,[3,[l is a widely- 
used framework in condensed-matter physics and quan- 
tum chemistry to calculate properties of many-electron 
systems, based on the simple local density instead of the 
less tractable A^-body wave function. 

One of the pillars of DFT is the Hohenberg-Kohn (HK) 
theorem Q , which in its original form proves that for any 
non-degenerate system of Fermions or Bosons [l| put into 
a local external potential v^^tir), there exists a unique 
functional of the local one-body density p(r) that gives 
the exact ground-state energy when p(r) corresponds to 
the exact ground-state density. A thorough mathemat- 
ical analysis of the foundations of the HK theorem was 
given by Lieb 0]. A crucial point is that as the theorem 
is based on the Ritz variational principle, it is valid only 
for systems described by a normalisable wave function [6| , 
i.e. for which a bound (ground) state exists. Various ex- 
tensions of the HK theorem have been proven, for exam- 
ple to spin-density energy functionals, for non-local exter- 
nal potentials, for relativistic, for time-dependent, or for 
superconducting systems [Sj]. The Kohn-Sham (KS) 
scheme furthermore provides a straightforward method 
to compute self-consistently the ground-state density in 
a quantum framework, defining the local single-particle 
potential (i.e. the non-interacting system) which repro- 
duces the exact ground-state density through an auxil- 
iary product state. 

The self-consistent mean-field (SCMF) approaches us- 
ing effective interaction that are widely used to describe 
the low-energy structure of atomic nuclei 0], resemble 
a KS scheme in many ways. Originally conceived as 
Hartree-Fock (HF) of Hartree-Fock-Bogoliubov (HFB) 
method based on an effective in-medium interaction, 
this framework has often been characterized as "nuclear 
DFT" [i, i, M, [m, E El, [3- The similarities be- 
come particularly obvious when the effective interaction 
is explicitely constructed as an energy functional depend- 
ing on various local densities and currents [1^1 • There 
are, however, important conceptual differences that pre- 
vent the straightforward mapping of the existing nuclear 
SCMF schemes onto the standard KS formalism for elec- 



tronic systems. For example, nuclei are self-bound, the 
intrinsic nuclear states obtained by SCMF methods of- 
ten break several symmetries of the nuclear Hamiltonian, 
and many extensions of the nuclear SCMF method aim 
at the explicit calculation of correlation effects instead of 
absorbing them into the functional. The present article 
addresses the first of these points, aiming at a KS scheme 
for self-bound systems. Similar efforts leading to approx- 
imate KS schemes have been made before 19, 25]. Here, 
we propose an alternative demonstration of the HK the- 
orem that carefully considers the technical issues arising 
from the separation of internal and center-of-mass co- 
ordinates required for for its application to self-bound 
systems, and that leads to an internal KS scheme.^ 



II. THE PROBLEM 
A. The role of the external potential 

In electronic systems, the wave function and density 
are defined in the frame attached to the center-of-mass 
(cm.) of the atomic nuclei. The latter also provide 
naturally the external potential Vextij), whose presence 
is compulsory to bind electrons that repel each other. 
The key point of the HK theorem is that the pure elec- 
tronic problem is universal, whatever the external field 
(provided it gives a bound state). 

In self-bound systems (as atomic nuclei or He droplets) 
the situation is intrinsically different because the net 
Fcrmion-Fermion (or Boson-Boson) interaction is attrac- 
tive. Thus, external fields are not necessary to obtain 
bound states, such that we are immediately in the cor- 
responding "pure" system, with the big difference that 
such systems physically exist. The absence of an ex- 
ternal potential, however, has as a consequence that the 



^ We follow here the nomenclature of [23j to refer to coordinates 
independent of the center-of-mass as "internal" ones, whereas 
we reserve " intrinsic" for symmetry-breaking states out of which 
bands of rotational states and/or parity vibrations can be mod- 
eled. 
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modeling of isolated self-bound systems is plagued by a 
center-of-mass (cm.) problem. For any stationary state 
with arbitrary total momentum P, the cm. will be delo- 
calized and evenly distributed in the whole space. Even 
more critical is that such laboratory wave function is not 
normalizable. This prevents any attempt to formulate 
DFT for isolated self-bound systems in terms of the lab- 
oratory density by simply taking the limit v^xt (r) — > in 
the HK theorem. Indeed, this density is an indeterminate 
constant [13, [111 , which forbids to construct a universal 
functional from it. It is of course the "internal density" 
Pint, i-G. the density relative to the system's cm., which 
is of interest. But standard DFT concepts as formulated 
so far are not applicable yet in terms of a well-defined 
internal density. 



B. The center-of-mass problem 

A second key point is that in a Hamiltonian and wave- 
function based description of an isolated self-bound sys- 
tem, the Hamiltonian should be explicitely translation- 
ally invariant to ensure Galilean invariance of the wave 
function ^. Hence, the iV-body wave function ip can be 
separated into a wave function F that depends on the 
position R = i J2f=i '^^ the cm. only, and an "inter- 
nal" wave function ipint that depends on the remaining 
{N — 1) Jacobi coordinates <^q, defined as = r2 — ri, 

6 = rg - . . . , ^N-l ^j^ivN- R) 

V'(ri, . . . ,rjv) - F(R) 7/.„t(6, • • • • (1) 

The r(R) describes the motion of the isolated system as 
a whole in any chosen inertial frame of reference (as the 
laboratory). The iptnt describes the internal properties 
and is function of the {N — 1) Jacobi coordinates. Of 
course it could also be written as a function of the TV 
coordinates , but one of them would be redundant [2^ . 

In this context pmt, rather than the laboratory den- 
sity p, thus becomes the natural quantity on which to 
construct DFT in a self-bound system. It is to be noted 
that for such a finite system it is impossible to construct 
a product state that has the required structure of 4'int- 
In a HF framework, one approximates instead directly 
il>(ri, . . . ,rpf) by a Slater determinant in N coordinates 
Tj in the cm. frame of the system. Consequently, the 
HF state contains (at least) one redundant coordinate, 
that introduces a spurious couplin g be tween the internal 
properties and the cm. motion [20|. For this reason. 



the HF approximation sacrifices "Galilean invariance for 
the sake of the Pauli principle" , to quote [llj . A rig- 
orous remedy is to perform projected HF, where projec- 
tion before variation on cm. momentum restores Galilean 
invariance at the price of abandoning the independent- 
particle description [2l|, [111- This reasoning does not 
hold, in principle, for DFT, where the key ingredient is 
the density, not an explicit iV-body wave function. 

A demonstration of a rigorous internal HK theorem has 
been made recently in [ij, (by two different ways), 
aiming at the correct separation of the internal proper- 
ties from the cm. motion, but none of them led to a 
rigorous internal KS scheme. A source term coupled to 
the TV-body internal density operator was introduced in 
[ill, allowing to express the exact total energy of a self- 
bound system as a functional of this operator. A scheme 
to construct a corresponding non-interating system in a 
systematic manner was proposed, but its link with the 
KS scheme of traditional DFT remains unclear. In [2^, 
it was shown that the internal energy of a self-bound sys- 
tem can be written as a functional of the internal one- 
body density and an approximate KS scheme was pro- 
posed, valid only if the cm. coordinate is treated as an 
adiabatic variable. A different approach to the problem 
is taken in Refs. [2l,[23|, where a (non-translationally in- 
variant) oscillator potential that traps the center-of-mass 
is added to the self-bound Hamiltonian. This aproach has 
the particular characteristic that it does not affect the 
internal properties of the system: the ground-state wave 
function is a wave packet that factorizes into the form of 
Eq. ([1]), r(R) now being a Gaussian and with that nor- 
malizable. The laboratory density p is then well defined 
and a KS scheme for p can be rigorously set up. The 
internal density pint can be deduced from p by deconvo- 
lution. However, the thus obtained energy functional 
and KS equations are neither an internal energy func- 
tional nor internal KS equations. Thus, the question of 
a rigorous formulation of an internal KS scheme compa- 
rable to SCMF calculations using an effective interaction 
remains open. Here, we propose a complementary way 
than those found in [H, [2^ to demonstrate the internal 
HK theorem, whose link with the traditional HK theo- 
rem is more clear. This directly leads to the formulation 
of a general internal KS scheme. 



III. DFT IN INTERNAL DEGREES OF 
FREEDOM 



^ Translational invariance, which states that the observables do 
not depend on the position of the cm., is a necessary, but not 
sufficient condition for the more fundamental Gahlean invari- 
ance, which ensures that observables are the same in all inertial 
frames. In case of a relativistic description of the quantum N- 
body system , Lorentz invariance has to be considered instead 
of Galilean invariance. 



A. Separation of internal and cm. coordinates 

We start from a general translationally invariant A^- 
body Hamiltonian composed of the usual kinetic energy 
term, and a two-body potential u which describes the 
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Fermion-Fcrmion (or Boson-Boson) interaction 



H 



N 2 

E Pi 
2m 

i=l 



N 



(2) 



i,j=i 
i>j 



For the sake of simplicity of the demonstration we as- 
sume a momentum-independent 2-body interaction and 
N identical particules. The generalization to 3-body in- 
teractions is straightforward, the generalization to sys- 
tems containing different types of particles will be dis- 
cused elsewehere. We rewrite the Hamiltonian using the 
iV — 1 Jacobi coordinates ^q, to decouple the internal 
properties from the cm. motion. The are to be 
distinguished from the N "laboratory coordinates" r^, 
and the "cm. frame coordinates" (r.^ — R) relative 
to the total cm. R. One can then separate ^ into 
H = HcM + H^nt, where Hcm = -ihy2M)AR (with 
M — Nm being the total mass) is a one-body operator 
acting in R space only, and Hint is a (N — 1) body op- 
erator in the ^q, space. Hint contains the interaction u 
which can be rewritten as a function of the [which we 
denote w({Cq}) for simplicity], and the internal kinetic 
energy, which is expressed in terms of the conjugate mo- 
mentum Tq, of and the corresponding reduced mass 
f^a = "^2^^- As [Hcm, Hint] = 0, the eigenstate of H 
can be build as a product of the form ((T|) with 



2M 

Hinti^int 



Eint'^int ; 



where 
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a— 1 



(3) 
(4) 



(5) 



There is no bound state for r(R) as the solutions are ar- 
bitrary stationary plane waves, leading to arbitrary cm. 
energy Ecm = /i^K^/(2M) and delocahzation of R. 
We will come back to the interpretation of r(R) below. 
By definition of a self-bound system, '\\}int is a bound, 
thus normalizable, state. The correponding total energy 
E = Ecm + Eint splits into 



E^'^\i,.M\ 

Eint \^ int J 



h^K^ ^ E,nt[lp,nt] 

2M irpint\i>int) 

(ipintlHintltpint) 



(6) 
(7) 



where the internal energy is obviously a functional of 
Tpint, and the cm. energy is parametrized by an arbi- 
trary K. We sec that the cm. properties (given by K) 
and the internal properties (given by ipint) are fully de- 
coupled. The ground state ipint of Hint is obtained by 
minimization of the total energy iJ^^^ [4'mt] for a given 
K, or equivalently of Eint['<pint] imposing normalisation. 

The previous steps have allowed to uniquely identify 
and separate the cm. motion. In traditional electronic 



DFT the problem does not show up as the electronic 
properties are defined in the frame attached to the cm. 
of the nuclei, where the nuclear background is accounted 
for by introducing an external local one-body potential 
Vext{j^) that provides the key ingredient of the HK the- 
orem. In the case of a self-bound system, v^xt is not 
compulsory. To faciliate the proof of the HK theorem, 
however, we introduce an arbitrary potential Vaux, which 
serves as a mathematical auxiliary and can be safely 
dropped at the end to recover an isolated self-bound sys- 
tem. 



B. A translational invariant auxiliary potential 

To conserve the separation of cm. and internal prop- 
erties, we cannot simply use a one-body potential of 
the form Ua«x(r). The potential Vaux should necessarily 
verify two conditions: (1) translational invariance and 
(2) as we are interested in the internal properties, the 
redundant cm. coordinate should be removed (as dis- 
cussed previously). These two conditions impose the 
form X^il^i ^aux (rj — R) as already used in [l^, [2^ , which 
corresponds to an arbitrary potential seen in the cm. 
frame. It can be expressed as a function of the Jacobi 

coordinates only, YliLi Wa«x(ri-R) = Vaux{{^a})] hence, 
it does not couple to cm. properties, the decomposition 
((H) for V still holds with Hint —>■ H,nt + Vauxii^a}) in 
(|4]). Of course, the associated internal wave function is 
modified accordingly, and consequently all internal ob- 
servables, but for the sake of simplicity we keep the same 
notations {ipint , Eint , Hint ) • 

As the next step, we evaluate the contribution of the 
auxihary potential term {ipint\vaux{{^a})\ipmt) to the in- 
ternal energy. First, we note that for any operator / that 
can be expressed through the Jacobi coordinates in po- 
sition representation (we note /({Ca}) when expressed 
through the Jacobi coordinates and /({r^}) when ex- 
pressed through the laboratory coordinates), we have the 
relation 



ilpint\fi{Ca})\lpint) 

Vi---d^iv-i^L({e4)/(Ua})V'»t({e4) 



dR dCl ■ • • d^N-l <5(R) V'*„t({ea})/({a}) V'mt({Ca}) 

= /dri...dr^<5(R)V*„t({r.})/({r,})V'™t({r.J). (8) 



We see that the "internal mean values" calculated with 
Ipint expressed as a function of the (A^ — 1) {^q} can 
also be calculated with ipint expressed as a function of 
the N coordinates {r,}. The transformation from the 
{£,a} to the {vi} introduces a <5(R) that represents the 
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dependence of the redundant coordinate on the others 

For the mean value of the auxiHary potential, relation 
(El) leads to 



N 



= I dri---dvN (5(R)|V'mt({ri})P ^^Vauxi^^t - R) 

i=l 

i=l 

X(5(77-(r, -R)) 



dr Vaux (r) Pint (r) 
where we have introduced the internal density 

P^ntir)/N 

= jdvi - - dvN (5(R)|V',nt({rJ)p 6{v ~ (r, - R)) 
= y"dri • • • drAT <5(R)|V™t({rJ)p <5(r - [vj, - R)) 
= [-^^^ ^d^di^---di^^^m) 



(9) 



(10) 



The density Pmt(r) is normalized to TV and a function 
of the cm. frame coordinates. The laboratory density 
is obtained by convolution with the cm. wave function 
as in [H, [i^l . The potential Vaux that is N body with 
respect to the laboratory coordinates (and (A^ — 1) body 
when expressed with Jacobi coordinates), becomes one 
body (and local) when expressed with the cm. frame 
coordinates. The energy Eint\i>int\ ©j and thus the total 
energy i?^^-'[V'mt] ®, are then to be complemented by 



i^,.-i?„.+/dr._(r)p.„.(r). 



(11) 



/dr Wa«x(r)pi„t(r), and as the ground state of ffi„i is 
obtained by minimization of Einti we can directly apply 
the usual proof of the HK theorem [3| and claim that 
for a non-degenerate ground state i/'mt, and for a given 
Fermion or Boson type (i.e. a given interaction u), the in- 
ternal energy Eint of a self-bound system, Eq. pip , can 
be expressed as a unique functional of pint- As already 
emphasized, the HK theorem is valid only for arbitrary 
"external" potentials that lead to bound ground states 
[|. As a direct consequence, the internal DFT scheme is 
valid only for potentials Vaux that lead to bound internal 
ground states -fAint- As for self-bound systems, described 
by our formalism at the limit —>■ 0, -tpint should by 
definition be a bound ground state, the previous conclu- 
sions still hold without breaking the consistency of the 
scheme. 



D. The internal Kohn-Sham scheme 

To recover the "internal" KS scheme, we assume, as 
in the traditional KS scheme, that there exists, in the 
cm. frame, a local single-particle potential (i.e. a A'^- 
body non-interacting system) which reproduces the den- 
sity pint of the interacting system. We develop pmt in the 
corresponding basis ifl^^^ of one-body orbitals expressed 
in cm. frame coordinates r ^ 



N 



Pint 



(r)=Ek™*«l 



(12) 



The KS assumption implies ^IntiPint] [l|; hence, we can 
rewrite Eint as ^ 



N 



Eint [P 



ExC [Pi: 



Y.(^lnt\^Mnt) + EH[Pint] (13) 
i=l 

+Exc[Pint] + Jdr 

'^aux ('^) Pint (l*) 

i j dv dr' jint (r, r') u{r - r') - Eh [pmt] 



N-l 9 



N 



C. The Hohenberg-Kohn theorem 

The internal energy Eint remains obviously a func- 
tional of ipint- As in its definition enters an arbi- 
trary one-body potential in the cm. frame of the form 



More generally, we can introduce a 5{R — a) where a is an ar- 
bitrary translation vector, which represents the position of the 
system's cm. in .the laboratory coordinates {ri}. For sake of 
simplicity of the notation, we chose a = without loss of gener- 
ality. 



* Even if only (A^ — 1) coordinates are sufficient to describe the 
internal properties, we still deal with a system of N particles. 
Thus, we have to introduce A'^ orbitals in the KS scheme if we 
want them to be interpreted (to first order) as single-particle 
orbitals and obtain a scheme comparable to SCMF calculations 
using effective interactions. 

^ To keep close contact with standard DFT, we make here 
the usual separation of the energy into direct (Hartree) and 
exchange-correlation parts. Owing to the complexity of the 
nucleon-nucleon interaction in the vacuum, strong correlations 
in the nuclear medium, and the appearance of three-body forces, 
it common practice in nuclear applications to construct approx- 
imate expressions for the entire functional. This, however, does 
not affect the conclusions of the present paper. 
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where we added and subtracted the internal Hartree en- 
eigy EniPint] = ^ Jdr dr' p,nt{r) Pint{r') u{r - r') for 
the direct part and the non-interacting kinetic energy 



'^iLii'filntl^f^l'Pint) ■ F'o'^ convenience, we introduced the 
local part of the two-body internal density matrix 



7mt(r,r') 

= Jdri---drN S{R)\ipint{{ri})f 

x^(r-(r, -R))<5(r'-(rj^, -R)) 
N{N -I) /N 



(14) 



N 



N 
N - 1 



d^i • • • d^N-3 

2 



nh ( f C r' + (jV-l)r Nr' \ 

Wintykli ■ ■ ■ An-3, — — , jj-rj j 



jV-1 f 



(using similar steps as in (fTO|) and — R 
and rAr_i - R = %E^£.n-2 - %^ )• Imt has the required 
normalisation to N{N — l)/2, is a function of the cm. 
frame coordinates and gives the two-body density matrix 
7(r,r') in the laboratory by convolution with the cm. 
wave function r(R). Applying Eq. ([8]) to the w({^q}) 
part of Hint and inserting Eq. gives directly 

{iJznt\u{{^a}Mint) = ^ [ dr dr' 7i„t(r,r') u{r - r') . 



(15) 

Following similar steps as in Eq. ([5]), one can show that 
the interacting kinetic energy can be rewritten as 



N-l ^2 

a—1 



(16) 



b2A N-l 2 



2M ^ 

Q — 1 



TV 



Jdr,--- dr^SiR) V*„t({r,J) ^ |^^™t({r,J) . 



Equation pB]) makes it clear that the dif- 
ference to the non-interacting kinetic energy 

EZiI drfi;A^)^^^l^ti^) comes, on the one hand, 
from the correlations neglected in the traditional 
independent-particle framework, but also from the cm. 
correlations (the (5(R) term in the previous expres- 
sion). The inclusion of the cm. correlations in the 
functional is the main difference to the traditional KS 
scheme. Still, the key point is that the internal pure 
exchange-correlation energy Exc[pint] can be expressed 
as a functional of pint- Varying Ej^[pint], Eq. or 
equivalently Eint[pint], Eq. ((T31), with respect to (/9-*j, 
and imposing normality of the {ipl„t}. 



(^Eint[pint] 



N 

E 



0, (17) 



leads to "internal" Kohn-Sham equations for the {ifllnt} 

- 2^^^ + UniPtnt] + UxciPznt] + Vauxjvlnt = ^iV\nt 

(18) 

with ?7H[pmt](r) = (5£'/i-[pi„i]/(5pi„t(r) and 
Uxc[Pint]ijc) = 5Exc[Pint]/5pint{i^), which is local as 
expected. Equation ([T51) has the same form as the tradi- 
tional KS equations formulated for non-translationally 
invariant Hamiltonians [iGj . Here, however, we have 
justified its use in the cm. frame for self-bound systems 
described with translational-invariant Hamiltonians, and 
shown that the functional form of Uxc[Pi7it\ differs by 
the inclusion of cm. correlations. Together with Eqs. (O 
and Eq. defines completely the total energy 

E-K.[pint\ as the sum of the cm. kinetic energy and of 
the internal energy. 



E. The laboratory density 

It is instructive to calculate the laboratory density p. 
Following Ref. [2^, one obtains 



p(r) = y"dR|r(R)|2p„,(r-R) 



(19) 



As r(R) is a plane wave, p(r) is constant. This confirms 
that the usual definition of the " laboratory density" lacks 
a meaningful interpretation for isolated self-bound sys- 
tems. Of course this full delocalisation does not occur in 
an experiment because observed self-bound systems are 
not isolated anymore. The observables related to r(R), 
i.e. position, momentum or kinetic energy, are used to 
transform all observables into the cm. frame, thereby 
explicitely using the Galilean invariancc. The key point 
is that the decoupling of the cm. motion allows to deduce 
internal properties preserving Galilei invariance, whereas 
r(R) is left to the choice of experimental conditions. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have shown in a way complementary 
to those proposed in Refs. [l^, [1^ that the total energy 
of a self-bound system can be expressed as a functional 
of the total one-body internal density pint ■ The energy 
in the laboratory frame contains the cm. wave vector K 
as a parameter that can be freely chosen according to ex- 
perimental conditions. Then, we have shown rigorously 
that the internal properties of the system are described 
by an internal KS scheme. The key difference to the tra- 
ditional HK/KS functional is the inclusion of cm. corre- 
lations. The question about the universal validity of the 
Kohn-Sham hypothesis, known as the "non-interacting 
v-representability" problem [l|, however, remains to be 
answered, as in traditional DFT. The internal DFT 
scheme proposed here provides a justification for the ap- 
plication of DFT to isolated ^Hc and '^He droplets 
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The present paper establishes also the first step towards 
an Kohn-Sham scheme applicable to nuclear structure 
physics. Further necessary developments are the gener- 
alization to two (or more) species of interacting particles, 
and the treatment broken rotational and space-inversion 
symmetry that requires to formulate the theory in term 
of the so-called " intrinsic" one-body density as defined in 
Ref. 29]. 
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